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NATIONAL ADVISORY COMMITTEE FOE AERONAUTICS 


RESEARCH MEMORANDUM 
for the 

Air Materiel Command, U* S. Air Force 


A THEORETICAL INVESTIGATION OF THE DYNAMIC LATERAL STABILITY 
CHARACTERISTICS OF THE MK -838 (EB-51) AIRPLANE 
By John W. Paulson 


SUMMARY 


At the request of the Air Materiel Command., TJ. S. Air Force, a 
theoretical study has teen made of the dynamic lateral stability charac- 
teristics of the MK -838 (XB-51) airplane. The calculations included the 
determination of the neutral-oscillatory-stability "boundary (R = 0 ) , the 
period and time to damp to one-half amplitude of the lateral oscillation, 
and the time to damp to one-half amplitude for the spiral mode. Factors 
varied in the investigation -were lift coefficient, wing incidence, wing 
loading, and altitude. 

The results of the investigation showed that the lateral oscillation 
of the airplane is unstable below a lift coefficient of 1.2 with flaps 
deflected 40° but is stable over the entire speed range with flaps 
deflected 20° or 0°. The results showed that satisfactory oscillatory 
stability can probably be obtained for all lift coefficients with the 
proper variation of flap deflection and wing incidence with airspeed. 
Reducing the positive wing incidence improved the oscillatory stability 
characteristics. The airplane is spirally unstable for most conditions 
but the instability is mild and the Air Force requirements are easily met. 


INTRODUCTION 


At the request of the Air Materiel Command, U* S. Air- Force, a 
theoretical study has been made of the dynamic lateral stability charac- 
teristics of the MK-838 (XB- 51 ) airplane. A three -view sketch of the 
airplane is shown in figure 1. The calculations were made at the Langley 
Laboratory on the Bell Telephone Laboratory * s 2-66744 relay computer. 

The analysis was made by the Langley free -flight -tunnel staff. 
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Calculations were made to determine the neutral-lateral-oscillatoiy- 
stability "boundary (E = 0), the period and time to damp to one -half amplitude 
of the lateral oscillation, and the time to damp to one -half amplitude for 
the spiral mode for several conditions of the airplane. The factors varied 
included lift coefficient, wing incidence, wing loading, and altitude. The 
results of the investigation are presented in the form of stability charts 
where the E = 0 "boundaries are plotted as functions of the effective - 
dihedral parameter -C and the directional-stability parameter C n p. 


The period and time to damp to one -half amplitude are presented as functions 
of lift coefficient for the various conditions investigated- 


SYMBOLS AND COEFFICIENTS 


S wing area, square feet 

c mean aerodynamic chord, feet 


V airspeed, feet, per second 

b wing span, feet' 

q dynamic pressure, pounds per square' foot 


p air density, slugs per cubic foot 

W weight, pounds 

I 

g acceleration of gravity, feet per second per second 

m mass, slugs (w/g) 

H-jn, relative -density factor based on wing span (m/pSb) 

i^. wing incidence, degrees 

a angle of attack of reference axis (fig. 2), degrees 

•q angle of attack of principal longitudinal axis of 

airplane, positive when principal axis is above 
flight path at the nose (fig* 2), degrees 

g angle between reference axis and principal axis, 

positive when reference axis is above principal 
axis at the nose (fig. 2), degrees 

angle between reference axis and horizontal axis, 
positive when reference axis is above horizontal 
axis at the nose (fig* 2), degrees 


Q 
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7 

t 

0 

4 

E 


k X 0 

k Zo 

K X 0 



% 




%Z 


C L 

Cn 

Cl 

°Y 

% 


angle of flight to horizontal axis, positive in a 
climb (fig* 2) , degrees 

angle of yaw, degrees or radians 

angle of sideslip, degrees or radians 

angle of hank, radians 

South's discriminant (E = BCD - AD 2 - B 2 E where A, 

B, C, D, and E are constants representing coefficients 
of the lateral- stability equation) 

radius of gyration about principal longitudinal axis, feet 

radius of gyration about principal vertical axis, feet 

nondimensional radius of gyration about principal 
longitudinal axis 

nondimensional radius of gyration about principal 
vertical axis ^k Zo /b^ 

nondimensional radius of gyration about longitudinal 
stability axis 


Co 2 cos 2 q + K Zo 2 sin 2 q J 
nondimensional radius of gyration about vertical stability 


axis 


N/^o 


P p P p 

cos T) + Kg sin q 


nondimensional product -of -inertia parameter 
((*Z 0 2 - 
lift coefficient (Lift/qS) 

yawing-moment coefficient (Yawing moment/qSb) 

rolling-moment coefficient (Soiling moment/qSb) 

lateral-force coefficient (Lateral force/qS) 

rate of change of lateral-force coefficient with 
angle of sideslip, per degree or per radian, as 

specified (^C y/^3) 

0mSSBSSffB9 


) 


cos q sin 
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rat© of charge of yawing-moment coefficient with angle 
of sideslip* per degree or per radian* as specified 

Q® n /Sp) 

rate of change of rolling-moment coefficient with angle 
of sideslip* per degree or per radian* as specified 

(Scj/se) 

rate of change of lateral-force coefficient with rolling- 
angular -velocity factor* per radian 



rate of change of rolling-moment coefficient with rolling- 


angular -velocity factor* per radian 



rate of change of yawing-moment coefficient with rolling- 
angular -velocity factor* per radian 

rate of change of rolling -moment coefficient with yawing- 
angular -velocity factor* per radian 

rate of change of yawing-mcment coefficient with yawing- 
angular -velocity factor* per radian 

rate of change of lateral-force coefficient with yawing- 
angular -velocity factor* per radian 

tail length (distance from center of gravity to rudder 
hinge line)* feet 

height of center of pressure of vertical tail above 
fuselage axis, feet 
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P 

r 


3>b 


s b 

X 


rolling angular velocity, radians per second 
yawing angular velocity, radians per second 
differential operator (d/ds-j^) 
distance along flight path, spans (vt/b) 
complex root of stability equation (c ± id) 


t 


time, seconds 


P 


period of oscillation, seconds 




time for amplitude of oscillation or spiral mode to 
change by factor of 2 (positive value indicates a 
decrease to half -amplitude, negative value indicates 
an increase to double amplitude) 

cycles for amplitude of oscillation to change by a 
factor of 2 


EQUATIONS OF MOTION 


The nondimensi onal lateral equations of motion (reference l) , 
referred to a stability-axes system (fig. 3 ) > are: 

Boll 


2 ^(% 2 »b 2 0 + KJCZV 1 *) “ + I + I CljV 

Ya w 

+ %zv¥) = CnpP + \ + \ 

Side slit 

2Hfc(V + Dbt) = + | Cy Vtf + CjJ + | Oj^DbV + (c L tan 7 ) + 

When ^ 0 e^ s ^> is substituted for $, i}r 0 e St) for \}r , and 3 Q e^ s b 
for 3 in the equations written in determinant form, X must be a root of 
the stability equation 

AX 1 *- + BX 3 + CX 2 + DX + E = 0 
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where 

A = . ]3 3 (kx 2 Kz 2 - Kxz 2 ) 

B = -2^- b 2 ^2K x 2 K z 2 Gy j3 + + K Z 2 C - ^xz, 2 ^ " %Z C Z r ' K XZ C n^) 

C = l^^C^Ojp + W^np + K Z 2c l v °Tp + 2 %% " K XZ C * r % 

" ^ifcPZp " C n p K XZ C Yp " | C n p C Z r + K XZ C n p C X p 
" K Z 2c y p G Zp " K X 2c Y r C np + K JZ G I r C l^ 

D = " ^ c n r c l p °rp _ ^b C Z p C np + | C np C 2r C Yp + ^b G n p C Zp 

+ 2[I 'b C L K XZ C np " 2|J ”b G L K Z 2c Zp " C np C L tan 7 

+ 2^^^ tan 7 + £ ’ £ %%% " t %%% 

+ t C “r C *J3% 

E - 1 ^(%% ■ %%) + 1 °L tan 7 (%% ' %%) 


The damping and period of the lateral oscillation are given by the 

equations T^yg ~ ~ P = ~^ and Ci/2 - where c and d are 

the real and Imaginary parts of the complex root of the stability equation. 
The damping of the spiral mode is determined similarly from one of the two 
real roots (usually the less stable one) of the stability equation. 

The conditions for neutral oscillatory stability as shown in reference 2 
are that the coefficients of the stability equation satisfy Routh's dis- 
criminant set equal to zero 

R = BCD -AD 2 - B 2 E = 0 


and that the coefficients B and D have the same sign, 
sign of the coefficient B is determined by the factors 


In general, the 




-C 


V 


and -C- 


which appear in the predominant terms of B- Thus, B is positive in the 
usual case if there is positive weathercock stability and positive damping 
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in roll. Hence, the coefficient D must he positive if R = 0 is a 
neutral-oscillatory- stability boundary . 


SCOPE OF TEE CALCULATIONS 


Calculations were made of the neutral-oscillatory-stability boundary, 
the pe. od and time to damp to one-half amplitude for the lateral oscilla- 
tion, and the time to damp to one -half amplitude for the spiral mode for 
several conditions of the airplane. The effect of lift coefficient was 
determined for flaps deflected and retracted for weights of 53>000 and 
32,017 pounds. The effect of varying wing incidence frcm 9° to 2° and 
-5* 3° "was also determined for a weight of 32,017 pounds with flaps deflected 
and the effect of changing altitude from sea level to 35>000 feet was 
determined for flaps up at a weight of 53 .>000 pounds. 


The aerodynamic and mass characteristics used in the calculation of 
the boundaries are presented in table I. The values of e n g)tail off 


and 


(%) 


tail off 


were obtained from force -test data furnished by the 


Glenn L. Martin Company. The tail-off values of C^ 


P 


C rLp> C l r ’ c n r > 


and Cy^ were estimated from data obtained in the Langley stability 


tunnel on a wing similar to that of the MK -838 airplane. The total 
derivatives including the contribution of the tail were estimated from 
the equations presented with table I. 


The roots of the stability equation were computed to determine the 
period and time to damp to one -half amplitude for each condition for which 
a boundary was calculated. Presented in table II are the aerodynamic and 
mass characteristics used in computing the roots of the stability equations. 
The values are essentially the same as those given in table I except that 
the mass characteristics are presented in a different form and the total 
derivatives are given. The values of C^ and C 7 used in these 

n 0 L $ 

calculations were obtained from force -test data furnished by the Glenn L. 
Martin Company. 

The variation of these values with flap deflection is shown in the 
following table: 


8 f 

(<leg) 

% 

(per deg) 

. ( per deg) 

4o 

0.00523 

0.00400 

20 

. 00490 

. 00290 

0 

. 00454 

. 00175 
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These values -were assumed to remain constant with lift coefficient for a 
given configuration of the airplane on the basis of the force-test results 
available on a model of this airplane. It would be expected, however, that 
the airplane would show an increase in at the higher lift coefficients. 

EESULTS AM) DISCUSSION 


The results of the calculations are presented in figures 4 to 1 6 - All 
results are presented in terms of lift coefficient and figure 14 is presented 
for convenience in interpreting the resuits in terms of airspeed. The test 
point in figures 4 to 9 represents the C n p and -C^ values of the airplane 

as shown by force-test results. 


Neutral-Oscillatory-Stability Boundaries 

Flaps down 40° . - The effect of lift coefficient is shown in figure 4(a) 
for the airplane with flaps down at sea level, iw = 9°* and a weight 
of 53,000 pounds. It is seen that the airplane is oscillatorily unstable 
below a lift coefficient of 1.2 and becomes stable at lift coefficients above 
this value. The results of figure 4(b) are for the same condition as 
figure 4(a) except that the weight is 32 , 017 pounds. Here again the air- 
plane is unstable below Cp = 1.2 and stable at hig h er values- 

The effect of wing incidence on the lateral stability of the airplane 
is shown in figure 5 for a lift coefficient of 0.4 with flaps down, light 
loading, and sea-level condition. Changing the wing incidence from 9° 
to 2° has a relatively small effect on the boundary although it is shifted 
downward far enough to make the airplane slightly stable. A further change 
in wing incidence to -5*3° (this change brought the inclination of the 
principal longitudinal axis of inertia with respect to the flight path 
to 0°) resulted in the airplane being very stable. The s ma ll shift in the 
boundary occasioned by the incidence change from 9° to 2° can be accounted 
for by the fact that with this change the principal axis of inertia was 
moved from -14*3° to -T*3°* Previous investigations (reference 3) have 
shown that the greatest effect of changes in the principal axis occurs at 
angles around 0°. Thus, when the incidence was changed so as to bring the 
principal axis from ~7*3 to 0°, a large downward shift in the E = 0 
boundary was obtained. 

Flaps down 20° . - The effect of lift coefficient on the E = 0 boundary 
for the airplane at sea level, half flap, i v = 5*5°* and 53*000 pounds is 
shown in figure 6(a). The boundaries for a weight of 32*017 pounds are 
given in figure 6(b) . From these results it is seen that the oscillation 
of the airplane is stable at Cp, = 0.4 and the stability increases with 
increasing lift coefficient. 
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Flans up . - The effect of lift coefficient on the oscillatory stability 
of the airplane with flaps retracted at sea level, i-^ = 2°, and a weight 
of 53 * 000 pounds is- shown in figure 7(a)* The airplane appears to have 
satisfactoiy stability even at a lift coefficient of 0.1 and increasing the 
lift coefficient resulted in a pronounced increase in stability. The 
results of the calculations made to determine the stability characteristics 
of the airplane with flaps up at 35*000 feet are presented in figure 7(h). 
These results show the same general effect of lift coefficient as for the 
sea-level condition. 

The results presented in figure 8 show the effect of a change in 
wing incidence from 2° to 0° for the flap-retracted condition at an 
altitude of 35*000 feet and a weight of 53*000 pounds* The 2° incidence 
change resulted in a fairly large downward shift in the boundary because 
of the fact that the change in the principal axis occurred near 0° (- 3 *4° 
to -1.4°) . 


Modified airplane . - A boundary was calculated for the airplane with 
flaps deflected 40° to determine the effect of several mass and geometric 
changes proposed by the Glenn L- Martin Company. The maximum wing incidence 
was assumed to be 7° instead of 9°* the tail length was increased 59 inches 
by moving the vertical tail to the rear of the fuselage, the negative 
dihedral was increased from ~3° to -6°, and the inclination of the principal 
longitudinal axis of inertia was assumed to be -4° instead of -3°* The 
calculations were made for a weight of 40,200 pounds and a full-flap 
placard speed of 165 miles per hour which corresponds to a lift coefficient 
of 1.05* The results of these calculations are presented in figure 9(a). 

In order to show- the effect of these changes a boundary is presented in 
figure 9(b) which was obtained by interpolating between the flap-down 
results of figures 4(a) and 4(b). The results of figure 9(a) show that 
the geometric and mass changes caused the airplane to have a small degree 
of oscillatory stability, whereas in the original condition (fig. 9(b)) it 
was slightly unstable. The modifications had very little effect on the 
boundary and the improved stability was due primarily to a change in C np 

and which shifted the airplane point with respect to the boundary. 


Period and Damping of the Lateral Oscillation 

The period of the lateral oscillation as a function of lift coefficient 
is presented in figure 10 for the various conditions investigated. It is 
seen that the period increases rather rapidly as the lift coefficient is 
increased from low to moderate values and increases at a lesser rate at the 
higher lift coefficients. For a given lift coefficient the period increases 
slightly as the flap deflection is increased and the period is greatest for 
the light loading. There is very little effect of altitude at low lift 
coefficients and at higher lift coefficients increasing altitude results in 
a slight decrease in the period. 
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The reciprocal of the time to damp to one -half amplitude is plotted 
as a function of lift coefficient in figure H. The larger values of the 

factor indicate the greater damping. It is seen from this figure 

* 1/2 

that the time to damp is negative helow Cj, = 1.2 -with flaps deflected 

40°, which indicates instability for the airplane; whereas, the damping 
is positive for the 20° and 0° flap conditions. The effect of changing 
wing incidence from 9° to 2° and “5*3° at Cl = 0.4 is shown for the 40° 
flap condition. These changes would make the airplane stable although 
the damping is still probably unsatisfactory with 2° wing incidence. The 
damping increases with lift coefficient for 40° and 20° flap conditions 
with the greatest change occurring with the heavy loading. With flaps 
retracted the damping decreases in going from Cl = 0.1 to Cl = 0*4 and 
then increases at lift coefficients beyond 0-4- The damping was relatively 
high, however, at all lift coefficients. The flap-retracted condition also 
shows that the effect of going from sea level to 35,000 feet resulted in 
a large decrease in the damping although the airplane was still quite 
stable . 


Perhaps a clearer picture of the results of the investigation is 
shown in figure 12 where the reciprocal of the cycles to damp to one -half 
amplitude is plotted as a function of lift coefficient. The cycles to 
damp are determined by dividing the time to damp by the period. If the 


value of 


1 

c l/2 


is greater than 0.5, it means that the lateral oscillation 


damps in two cycles or less. From figure 12 it is seen that with flaps 
up, all conditions damp in two cycles or less; whereas with flaps down 
40°, only the heavy loading at the highest lift coefficients damps in 
two cycles. The low-lift -coefficient conditions with flaps up and an 
altitude of 35,000 feet are actually impossible to attain with this air- 
plane because of the high speed required as shown in figure 13 * 


From the results presented it is seen that there is a range of lift 
coefficients for oscillatory stability for the airplane for all flap 
conditions even though the range becomes quite small with a flap deflection 
of 40°. If care is taken to insure that the airspeed is low enough before 
a given flap deflection is reached, then the dynamic lateral stability 
might be satisfactory over the entire speed range. One suggested means 
for insuring the proper flap setting with airspeed has been an automatic 
flap -unloading device. This device would allow flap deflection and wing 
Incidence to change with airspeed such that a stable condition was always 
maintained. Figure 14 has been prepared to show how such a flap-unloading 
device might be applied. The reciprocal of the cycles to damp to one-half 
amplitude has been plotted as a function of airspeed in miles per hour for 
flap deflections of 0°, 10°, 20°, 30°, and 40°. (Th.e curves for 10° 
and 30° flap deflections were obtained by interpolation. ) A linear 
variation of flap deflection with airspeed was assumed and cross-plotted 
as dashed curves in figure 14 to show the damping of the airplane for 
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each flap position. The lower of the two curves is based on the assumption 
that the flap is released at 250 miles per hour and reaches full deflection 
at 165 miles per hour whereas the upper curve was "based on a flap release 
speed of 220 miles per hour and a full-flap speed of 155 miles per hour. It 
is seen that the airplane remains stable for all conditions although the 
oscillation requires more than two cycles to damp with 40° flap deflection 
at a speed of 165 miles per hour. This damping might not be objectionable 
in the landing condition and a reduction in speed to 150 miles per hour 
would bring the oscillation to a two-cycle damping. 

The fact that the lateral oscillation damps in two cycles does not 
necessarily indicate that the lateral stability is satisfactory. Some 
present-day airplanes when flying at high speeds have encountered undesirable 
short-period oscillations that damp in less than two cycles. The Air Force 
and Navy have therefore set up requirements which call for greater damping 
of the shorter -period oscillations. Presented in figure 15 is a plot showing 
how the MX- 83 8 airplane damping compares with the requirements of both the 
Air Force and NACA. (See references k and 5* ) It is seen that the Air 
Force requirement is much more stringent than the NACA requirement for 
periods from 1 to 3 seconds. Figure 15 shows that several conditions of 
the airplane that meet the NACA requirement would actually be unaatisfactory 
or very close to being marginal when considering the Air Force requirement. 


Spiral Stability 

The reciprocal of the time to damp to one-half amplitude for the 
spiral motion is presented in figure 1 6 for al 1 conditions investigated. 

It is seen that most conditions are spirally unstable but the instability 
is not severe. In all cases the spiral instability increases up to a 
lift coefficient of about 0.8 and then decreases at higher lift coefficients. 

Although neither the Air Force nor the NACA flying-qualities 
requirements call for spiral stability, they do require that any spiral 
instability present should be mild. The Air Force requirement is for the 
spiral motion not to double amplitude in less than 4 seconds. It is seen 
from figure 1 6 that the spiral instability of the MX-838 airplane will be 
acceptable for all conditions. 


CONCLUSIONS 


The following conclusions are drawn from the results of the theoretical 
study of the dynamic lateral stability characteristics of the MX-838 air- 
plane: 

1* The lateral oscillation of the airplane is unstable below a lift 
coefficient of about 1*2 with the flaps deflected ^0° but is stable over 
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the entire speed range -with flap deflections of 20° or 0°. Satisfactory- 
oscillatory stability can probably be obtained for all lift coefficients 
with the proper variation of flap deflection and wing incidence with 
airspeed. 

2. Beducing the wing incidence from 9° to 2° with flaps deflected 
40° makes the lateral oscillation of the airplane slightly stable. A 
further change in wing incidence to -5. 3° (this change brought the inclina- 
tion of the principal axis to 0°) results in the airplane being very stable. 

3. The airplane is spirally unstable for most conditions but the 
instability is mild and the Air Force requirements are easily met. 
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TABLE I 

FACTORS USED HI CALCULATING THE OSCILLATORY STABILITY BOUNDARIES 


Flap 

deflection 

(deg) 

B 

a 

(deg) 

n 

(deg) 

Q 

W 

(lb) 

( slugs /cu ft) 


B 

k 2 
"0 

(^p) tail off 

(^ n p) tail off 

( C *p) tail off 

( Cn p)tall off 

( Cl r)tail off 

( C «r)tail off 




-11.3 

-14.3 


53,000 

O.OO2378 

23.6 

5-37 

16.5 

-0.57 

-0.074 

-0.35 

-0.05 

0.22 

-0.125 



Hi 

-6.1 

-9.1 

9 

53,000 

.002378 

23.6 

5-37 

16.5 

-•57 

-.074 

-•35 

-.05 

.32 

-•15 



1.2 

-•9 

-3-9 

9 

53,000 

.002378 

23.6 

5-37 

16.5 

-57 

-.074 

-•35 

-•05 

•35 

-.175 



1.6 

4.3 

1.3 

9 

53,000 

.002378 

23.6 

5-37 

16.5 

-57 

-.074 

-•35 

-.05 

•35 

-.20 



.4 

-11.3 

-14.3 

9 

32,017 

.002378 

14.27 

6-56 

20.0 

-57 

-.074 

-•35 

-.05 

.22 

-.125 

40 


.8 

-6.1 

-9.1 

9 

32,017 

.002378 

14.27 

6*56 

20.0 

-57 

-.074 

-•35 

-.05 

•32 

-.15 



1.2 

-•9 

-3-9 

9 

32,017 

.002378 

14.27 

6.56 

20.0 

-57 

-.074 

-•35 

-.05 

•35 

-175 



1.6 

4-3 

1-3 

9 

32,017 

.002378 

14.27 

6. 56 

20.0 

-•57 

-.074 

-•35' 

-05 

•35 

-.20 



.4 

-4.3 

■-7-3 

2 

32,017 

.002378 

14.27 

6. 56 

20.0 

-•57 

-.074 

-35 

-.05 

.22 

-.125 



.4 

3-0 

0 

-5-3 

32,017 

.002378 

14.27 

6-56 

20.0 

-•57 

-.074 

-•35 

-.05 

.22 

-.125 



1 .05 

-.8 

-4.8 

7 

40,200 

.002378 

17-9 

6.02 

18.6 

-.57 

-.074 

-•35 

-.05 

•34 

-.165 



r.h 

-3-8 

-6.8 

5-5 

53,000 

.002378 

23.6 

5-37 

16.5 

-.45 

-11 

-•35 

-.05 

.18 

-.075 



.8 

1.5 

-1.5 

5-5 

53,000 

.002378 

23.6 

5-37 

16.5 

-.45 

-.11 

-•35 

-.05 

.28 

-•095 



1.2 

6-9 

3-9 

5-5 

53,000 

.002378 

23.6 

5-37 

16.5 

-.45 

-.11 

-•35 

-.05 

•30 

-.11 

20 


.4 

-3.8 

-•4.8 

5-5 

32,017 

.002378 

14.27 

6> 56 

20-0 

-.45 

-.11 

-•35 

-.05 

.18 

-•075 



.8 

1.5 

-1.5 

5-5 

32,017 

.002378 

14.27 

6*56 

20.0 

-45 

-.11 

-35 

-.05 

.28 

-.095 



J.2 

6.9 

3*9 

5-5 

32,017 

.002378 

14.27 

6.56 

20.0 

-.45 

-.11 

-•35 

-.05 

•30 

-.11 



f.l 

-.4 

-3-4 

2 

53,000 

.002378 

23.6 

5.37 

16.5 

-.34 

-15 

-•35 

-.05 

0 

-.015 



.4 

3-8 

.8 

2 

53,000 

.002378 

23.6 

5-37 

16.5 

-•34 

-15 

-•35 

-.05 

• 15 

-.025 



.8 

11 

6.2 

2 

53,000 

.002378 

23.6 

5-37 

16.5 

-34 

-.15 

-•35 

-.05 

.24 

-•035 



1.2 

B 1 

12.0 

2 

53,000 

.002378 

23.6 

5-37 

16.5 

-34 

-15 

-•35 

-05 

.25 

-.05 

0 


.1 

Hi 

-3-4 

2 

53,000 

.000736 

76-3 

5-37 

16.5 

-.34 

-15 

-35 

-.05 

0 

-.015 



.4 

3.8 

.8 

2 

53,000 

.000736 

76-3 

5.37 

16.5 

-.34 

-15 

-35 

-.05 

• 15 

-.025 



.8 

K3 

6.2 

2 

53,000 

.000736 

76-3 

5-37 

16-5 

-34 

-15 

-•35 

-.05 

.24 

-•035 



1.2 

IWg 

12.0 

2 

53,000 

.000736 

76.3 

5-37 

16.5 

-.34 

-.15 

-•35 

-.05 

.25 

-.05 



L-i 

m 

-1.4 

0 

53,000 

.000736 

76.3 

5-37 

16.5 

-.34 

-15 

-•35 

-.05 

0 

-.015 


Constants: 
b = 53-3 ft 
S = 550 sq ft 
7 s 00 

Z = 0.58 ft (0.67 

D 

y P = o 
Variables: 

(S)t.u ■ °- s “ - 1 - 0 

| 0 0.08 to 0.185 


Total 


( 0, r)% 


derivatives Including 
-0-35 

\u 

+ 2 


tall contributions 
off 



determined from the following equations: 

2 . 

7 « 00' " v -VVnfO 2 (b " b 8ln “} 7 (■‘Stall off (°*p)total = ( C } V )c n?a0 '*(*'* Bln *) 

N V 0 “ ( C ^r)tail off + 2 (b " b ( ' C np)tai 1 off (%)total “ ( C V)c^ 0 + 2 (b ‘ b 0ln «)% 

.0.05 + a(| - I sin a) (-C np ) taU ^ (%) total “ (S)c ap -o + S (f ' 7 ° ln *)% 

^b(‘%)tail off ( n f) total B ( Cn r)c n ^= 0 ' 2 (b %) 


l % 



& 


c 


0T381S °OM wa VQVM 






























TABLE II 


table nr 

FACTORS USED IN CALCULATING BOOTS OF LATERAL- STABILITY EQUATIONS FOB DETERMINATION OF PERIOD AND TIME TO DAMP 


£ 


Flap 

deflection 

(deg) 

C L 

4 

(deg) 



C?)‘ 

-% 2 

K Z 2 

K xz 

(^33) total 

(^ 7 ^) total 

(^r^total 

( Cn p)total 

(^ n r) total 

m 

pSV 



f 6.4 

- 14.3 

23.6 

0.0102 

0.0958 

0.0154 

0.0906 

-0.0180 

- 0.4466 

- 1.15 

0. 4247 

0.1547 

- 0.5588 

2.8 



.8 

- 9-1 

23.6 

• 0102 

.0958 

.0123 

•0937 

-.0127 

-.4134 

- 1.15 

.4858 

.1158 

-.5838 

3-94 



1.2 

- 3-9 

23.6 

• 0102 

.0958 

.0106 

.0954 

-.0058 

-.3869 

- 1.15 

.4765 

•0765 

-.6088' 

4.85 



1.6 

1-3 

23.6 

.0102 

.0958 

• 0102 

•0957 

.0019 

-•3675 

- 1.15 

.4371 

•0371 

-.6338 

5.60 



.4 

- 14.3 

14.27 

.0151 

• l 4 o 8 

.0228 

•1331 

-.0264 

-4466 

- 1.15 

.4247 

•1547 

-.5588 

2.17 


< 

.8 

- 9-1 

14.27 

.0151 

.1408 

.0182 

•1377 

-.0187 

-.4134 

- 1.15 

.4858 

.1158 

-.5838 

3-06 

40 


1.2 

- 3*9 

14.27 

.0151 

.1408 

.0156 

.1403 

-.0085 

-.3869 

- 1.15 

*4765 

• 0765 

-.6088 

3-7 6 



1.6 

1.3 

14.27 

.0151 

.1408 

.0152 

.1407 

.0028 

-•3675 

- 1.15 

.4371 

•0371 

-.6338 

4-35 



.4 

- 7.3 

14.27 

• 0151 

.1408 

.0172 

.1387 

-•0153 

-.4035 

- 1.15 

.3170 

.1022 

-.5588 

2.17 



.4 

0 

14.27 

.0151 

.1408 

.0151 

• l4o8 

0 

-•3717 

- 1.15 

•3346 

.0470 

-.5588 

2.17 



[l -05 

- 4.8 

17.9 

.0128 

.1218 

.0136 

.1210 

-.0090 

-.3869 

- 1.15 

.48a 

.096 

-.744 

3-94 



r .4 

- 6.8 

23.6 

• 0102 

.0958 

.0114 

.0946 

-.0098 

-.4031 

- 1.10 

• 3346 

.1046 

-.5275 

2.8 



.8 

- 1.5 

23.6 

.0102 

.0958 

.0103 

.0947 

-.0022 

-.3782 

- 1.10 

.3928 

.0628 

-.5475 

3-94 



1.2 

3-9 

23.6 

.0102 

.0958 

.0106 

• 0954 

.0058 

-.3610 

-1.10 

,3705 

.0205 

-.5625 

4.85 

20 

< 

.4 

- 6.8 

14.27 

.0151 

.1408 

.0169 

.1390 

-.0144 

-.4031 

-1.10 

•3346 

.1046 

-.5275 

2.17 



.8 

- 1.5 

14.27 

.0151 

.1408 

.0152 

.1407 

-.0033 

-•3782 

-1.10 

.3928 

. 0628 

-.5475 

3-06 



1.2 

3-9 

14.27 

.0151 

.1408 

.0156 

.1403 

.0085 

-.3610 

- 1.10 

•3705 

.0205 

-.5625 

3-76 



r.i 

- 3-4 

23.6 

•0102 

.0958 

• 0105 

.0956 

-.0050 

-.3881 

-1.05 

.1345 

.0845 

-.4906 

1*395 



.4 

.8 

23.6 

.0102 

.0958 

.0102 

.0959 

.0012 

-,3708 

-1.05 

.2496 

.0496 

-.5006 

2.8o 



.8 

6.2 

23.6 

.0102 

.0958 

.0112 

.0949 

.0090 

-.3564 

-1.05 

.2952 

.0052 

-.5106 

3-94 



1.2 

12.0 

23.6 

.0102 

.0958 

♦0139 

.0921 

•0159 

-.3502 

-1.05 

.2581 

-.0419 

-.5256 

4.85 

0 

<j 

.1 

- 3*4 

76-3 

.0102 

.0958 

.0105 

.0956 

-.0050 

-.3881 

-1.05 

.1345 

.0845 

-.4906 

2.51 



.4 

.8 

76.3 

.0102 

.0958 

• 0102 

•0959 

.0012 

-.3708 

-1.05 

.2496 

.0496 

-.5006 

5-02 



.8 

6.2 

76.3 

.0102 

.0958 

.0112 

.0949 

.0090 

-•3564 

-1.05 

• 2952 

.0052 

-.5106 

7-10 



1.2 

12.0 

76.3 

.0102 

.0958 

•0139 

.0921 

.0159 

-.3502 

-1.05 

.2581 

-.0419 

-.5256 

8.69 



.•1 

- 1.4 

76.3 

.0102 

.0958 

.0103 

.0958 

-.0021 

-.3792 

-1.05 

• 1179 

.0679 

-.4906 

2.51 



Cy r = 0 
7 = 0° 
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Figure 1 .- Three-view drawing of the MX-838 (XB-51) airplane 



Figure 2.— System of axes and angular relationship in flight. Arrows 
indicate positive direction of angles, t) =» Q — y — e. 
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Figure 3 •— The stability system of axes. Arrows indicate positive 
directions of moment s, forces, and control— surface deflections. 

This system of axes is defined as an orthogonal system having their 
origin at the center of gravity and in which the Z-axis iB in the 
plane of symmetry and perpendicular to the relative wind, the 
X-axis is in the plane of symmetry and perpendicular to the Z-axis, 
and the Y— axis is perpendicular to the plane of symmetry. 
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0 WZ . 004 .006 MS . 010 .otz 

Effective - dihedral parameter , - 


(a) 53*000 pounds. 

Figure 4.- Effect of lift coefficient on the E = 0 boundary of the MX-838 airplane. 

Flaps 1*0°; sea level. 
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Directional - stability parameter , 



0 .003 .004 . 006 .008 . 0/0 . 013 . 0/4 

Effective- dihedral parameter , - C\^ 


m 


(l) 32,017 pounds. 
Figure k,— Concluded. 
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Di reef. zonal' <sfahiJ j fy para/r?e ter 



0 .002, .004 . 006 ' 006 * 010 . 0/2 . 014 

Effective - dihedral parameter ; - Ca 


Figure 5 .- Effect of wing incidence on the E = 0 boundary of the MX-838 airplane. Flaps 40°; 

sea level; 32,017 pounds j Cj, = 0A. 
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Directional- stability parameter , 
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Effective - dihedral parameter ; - 

(a) 53,000 pounds. 

Figure 6.- Effect of lift coefficient on the E = 0 boundary of the MX -838 airplane. 

Flaps 20 ° } sea l evel. 
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(b) 32,017 pounds. 

Pigi^^^Oonoludsd . 
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D/recf/o/iaJ- sfabif/fy parameter 



Figure 7»~ Effect of lift 


of the MX -838 airplane. 
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Figure 8 .— Effect of wing incidence on the R = 0 "boundary of the MX — 838 airplane. 

35,000 feet; 53*000 pounds. 
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Figure 9,- Effect of modifications to the MX-838 airplane on the K = 0 boundary. 

sea level; 40,200 pounds; Cl = 1.05. 



Flaps 40°; 
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Figure 10.— Period of 













*••• 



Figure 12.- Damping of the lateral oscillation of the MX-838 airplane in 
terms of the reciprocal of the cycles to damp to half amplitude. 
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Figure 13.- Effect of altitude and weight qji the airspeed required for 

the MX -838 airplane. 


a given lift coefficient 
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Figure 14 .- Design chart shewing the effect of flap setting and airspeed on the damping of the 

lateral oscillation of the MX-838 airplane. 
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Figure 15 .— Comparison of the damping of the MX -838 airplane with respect 
to the HACA and Air Force requirements. 53,000 pounds. 
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Figure 1 6 .- Damping of the spiral motion of the MI -838 airplane in terms 
of the reciprocal of the time to damp to one— half amplitude. 






